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Abstract 

We study the production of C = + charmonium states X in e~''e~ — > 7 + X at B factories with 
X = r]c{nS) (n=l,2,3), Xc,j{fnP) (m=l,2), and ^^02(10). In the S and P wave case, contributions 
of QED with one-loop QCD corrections are calculated within the framework of nonrelativistic 
QCD(NRQCD) and in the D-wave case only the QED contribution is considered. We find that 
in most cases the one-loop QCD corrections are negative and moderate, in contrast to the case of 
double charmonium production e'^e^ ^ J /ip + X, where one-loop QCD corrections are positive 
and large in most cases. We also find that the production cross sections of some of these states in 
e~^e~ — > 7 + X are larger than that in e~^e~ — > J/i/; + X by an order of magnitude even after 
the negative one-loop QCD corrections are included. We then argue that search for the X(3872), 
X(3940), Y(3940), and X(4160) me~^e~ — > 7 + X at B factories may be helpful to clarify the nature 
of these states. For completeness, the production of bottomonium states in e~^e~ annihilation is 
also discussed. 
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I. INTRODUCTION 



In recent years there have been a number of exciting discoveries of new hidden charm 
states, i.e. the so called XYZ mesons, by Belle, BaBar, CLEO, CDF, and DO collabora- 
tions (for recent experimental and theoretical reviews and related references see Ref. 
Among the XYZ states, the charge parity C=+ states such as X(3872), X(3940), Y(3940), 
Z(3930), and X(4160) are particularly interesting and the interpretations for their nature 
are still very inconclusive (except for the Z(3930), which is assigned as the Xc2(2P) meson). 
The experimental results for these C=+ states have induced renewed theoretical interest 
in understanding the mass spectrum, decayandproduction mechanisms of charmonium or 
charmoniumlike states (see, e.g., Refs. . Among others, the double charmonium 

production in e^e~ annihilation at B-factoriesp, [7| turned out to be a good way to find the 
C=+ charmonium or charmoniumlike states, recoiling against the easily reconstructed 1 
charmonium J/ip and ^'(25'). In addition to the r]c,r]c{2S) and Xco, the X(3940) (decaying 
into DD*) and X(4160) (decaying into D*D*) have also been observed in double charmo- 
nium production. Since the quantum number of the photon is the same as J/ip, it will be 
interesting to see whether the C=+ charmonium or charmoniumlike states can be found in 
the process e+e^ ~^ 7* ~^ 7 + -^, where X is a C=+ state recoiling against thephoton. The 
production rates of such processes have been calculated at tree level in QED|8|. 

It has been known for some time that the one-loop QCD radiative corrections are very 
important in double charmonium production in e~^e~ annihilation. The observed double 
charmonium production cross section for e~^e J/ipVc is larger than the leading- 
order (LO) calculations in NRQCD^ by an order of magnitude flO|. and later it was found 
that these discrepancies could be largely resolved by the next-to-leading-order (NLO) QCD 
12| combined with relativistic corrections jisl. Therefore, it is necessary 



corrections 
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to examine whether the one-loop QCD (i.e., 0(0;^)) corrections are also important for the 
processes e~^e~ 7* ^ 7 + X. In fact, the one-loop QCD radiative correction to e~^e~ — *■ 
7* — s> 7 + ?7c has been investigated elsewhere |l5l. Il6 |. 

Another interesting point is about 1"^^ charmonium. At B factories the observed produc- 
tion cross sections in e+e~ annihilation to J/i'Tjc, ilj{2S)rjc^ J / iprjci^S) , J /ipXco, and ip{2S)xco 
are large, but no signals for J /ipXci,2 have been seen. This is in line with the calculations in 
NRQCDpil], in which the predicted production rates of J I^Xc\,2 are relatively suppressed. 
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We wonder whether the cross section of 1"'""'" charmonium (including Xci and its radial exci- 
tations) associated with a photon could be large in e~^e^ ^ 7* ^ 7 + X. If this is the case, 
we might have a chance to search for the Xci as well as the X(3872) in e~^e~ —^7* — > 7 + X, 
since the X(3872) could be a Xcii'^P) dominated state but mixed with some D^D*^ compo- 
nent, in one of the possible interpretations. This is also useful to the search for the Y(3940), 
which has been seen in the decay B Y{394:0)K followed by y(3940) J/ipu, and is also 
a possible candidate for the XciC^P) (or Xcoi'^P))- Of course, these states could have some 
more exotic nature, being molecules, tetraquarks, or charmonium hybrids. 

In this paper, we compute the QED (at tree level) and one- loop QCD (0(as)) corrections 
to the processes e^e~ — 7* — > 7 + X, where X are Xcj,Vc, and their radially excited 
states, all with charge-parity C = +1. We find the cross sections for rjc, its radial excited 
states and xa, Xci(2-P) are relatively large. Despite of the large background from initial 
state radiation (ISR), we still expect they could be seen in the 7 recoil spectrum with 
higher statistics in the future. The remainder of the paper is organized as follows. In Sec. 2 
we outline the QED calculation and some basic techniques for numerically computing the 
one-loop QCD correction. The QED and one-loop QCD corrections to cross sections for 
e^e^ ^ 7 + X at B factories are given in Sec. 4, and we also analyze and discuss our 
results. In the Appendix, we show some basic integration expressions. 

II. QED CALCULATION 




FIG. 1: The tree QED diagram for e~^e — > 7 + X 

The Feynman diagram for the exclusive process e^e^ — > 7 + X at order a^a^ is shown in 
Fig.l, where X is a heavy quarkonium with charge-parity C = +1, and there is another quark 
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line-flipped one. In the nonrelativistic limit, the factorization formula for heavy quarkonium 
production in the NRQCD framework is equivalent to that in the color-singlet model. And 
in our case, the amplitude A4 for e+e^ — j> 7 + X can be expressed as 

M(e-e- - 0^ + X) = E E E / (4%5^<«° - 

{LL,, SS,\JJ,){z,j\l)A{e+e- ^ 7 + Q^^(| + q) + QI^{^ - q)) (1) 

where P is the momentum of X state, 2q is relative momentum between Q and Q in the 
rest frame of X state, and {LL^; SSz\JJz), {si, S2\SSz) and = Sij/\^N^ are the spin- 

SU(2), angular momentum C-G coefficients and color-SU(3) C-G coefficients for QQ pairs 
projecting onto appropriate bound states respectively. And A is the standard Feynman 
amplitude denoting e+e" — » 7 + Ql^{^ + q) + Q^ii: ~ 

The Feynman amplitude part can be evaluated by introducing the spin projection 

operator 

p p 

PssXP^ q) = ^{su S2\SSz)v{— - g; S2)u{— + q; si). (2) 
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Expanding the operator in terms of the relative momentum q, we get the leading-order 
nonvanishing terms for the S-, P- and D-wave case respectively. The results of the spin- 
triplet and spin-singlet projection operators and their derivatives with respective to the 



relative momentum q^ are given below 



19|: 



PisAP.O) = ^r{Sz){^ + 2mQ), (3) 

Pi%iP,0) = + - (^-2mQ)f (^.)7l. (4) 

Poo(P,0) = ^75(/' + 2mQ), (5) 

Pp"q^(P,0) = ^ il^i^ - 2mQ)'y'^ + 7^(/' - 2mQ)7")7^ + vanishing terms (6) 
8v2mQ 

After integrating g", we get the amplitudes for S, P and Z^-wave heavy quarkonium produc- 
tion respectively: 

-M(7 + r]c) = I ^^oo(q)Tr[PooO]|,=o (7) 

A^(7 + Xcj) = J ^g«^i^(q)(LL,; JJ.)4(5.)Tr[Pf5^0" + P,%0]U=o (8) 



where O is the 7 matrix relevant to the Feynman amphtude A, and and O"^ are the 
first and second derivatives of O with respect to Qa respectively. 

The integrals of the wave function in momentum space are related to the radial wave 
function Rs{0), -Rp(O) and Rd{0) in coordinator space at the origin for the S—, P— and 
D— wave cases respectively: 

"^'"^ ^oo(q) = -^Rs{0), (10) 



I ^gaV'in^(q) = -^eUL,)^R'p{0), (11) 

where €°'{Lz) is the polarization vector of L = 1 (P-wave) system and £^^{Lz) is the polar- 
ization tensor of L = 2 (D-wave) system. 

For spin-triplet P-wave states, the projection of the L — S coupling of the spin vector 
€*{Sz) and orbital vector e*{Lz) onto total angular momentum J for J = 0, 1, 2 are 

e*(5,)e;(L,)(l,L,; 15,|0,0) = -^H,^, (13a) 

e;(5,)e;(L,)(l, L,; 15,|1, J,) = -^e"/^''-p,e*( J,), (13b) 
el{S,)e;{L,){l,W, 15,|2, J,) =£;^, (13c) 

where Ua/s — {—gap + ^r-)- For the total angular momentum J —1 and J = 2 states, the 
sums over all possible polarizations are given by 

eaiJzVpiJz) = n„/3, (14a) 

y^,^al3{Jz)£*a'0'{Jz) = l/2(nQ,Q,/n/j/3/ ITa/j/ riaZ/j ) - -UapUa'P'- (14b) 

With the help of the formula introduced above, we get the final QED analytic expressions 
for the exclusive process e+e~ — > + X: 

a(e-e- - 7 + ^c) = '"'^''"^f ^''^^ " / dn{l + cos\e)) (15a) 
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cr e e 



7 + Xco) 



3a^et\R'piO)\^{l-3rf 
s^m^(l — r) 



dn{l + cos\9)) 



(15b) 



cr e e 



7 + Xci) 



18a^ei\R'piO)\' 



dn{l + 2r + (1 - 2r) cos\9)) 



(15c) 



a e^e 



7 + Xc2) 



6a^et\R'piO)\' 
s^m^(l — r) 



y cil](l + 6r + 6r^ + (1 - 6r + Gr^) cos2(^)) (15d) 



a e e 



7 + '^2 



15a3e4|^^ 



dn{l + coi?{d)) 



(15e) 



where r = M^js^ cos{6) is the angle between J/?/' and the initial beam axis. For the cc 



system we set Mx = 2mc. If we replace 
results of ^Pj are consistent with those in Ref. 



by 



1 (Osepj)) 



jy2_i — 2J+1 — ' QED 
201. For bb states, the result can be obtained 



by changing Cc to e?,, rric to rrify and the values of the wave-functions for charmonium states 
to those for bottomonium states. 



III. ONE-LOOP QCD CALCULATION 

Now we proceed to calculate the one-loop QCD corrections. The numerical calculation of 
one-loop QCD corrections is performed with the help of Feyncalc and Looptools. At the 
one-loop level of QCD, there are eight Feynman diagrams. We show four of them in Fig.2, 
and the other four can be obtained by reversing the direction of the charm quark line. At 
order a^ag, the cross section for e~^e~ — > 7 + X is 

dcr OC \Mtree + MqcdI'^ 

= \Mtree\^+2ReiMl,,MQCD)+Oia^a^,), (16) 

where A4qcd means the one-loop QCD amplitude. The on-shell scheme is adopted and 
then the self-energy renormalization constant Zi and vertex renormalization constant Z2 
are chosen to be 

SZ^' = - J- + 7^ - 4 - A _ log(^), (17) 

£uv £iR 
5Z^^ = 5Z^^ (18) 
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FIG. 2: The one-loop QCD diagrams for e~^e 



7 + X 



where we omit the coefficient before the self-energy renormalization constant and part of 
the infrared divergence term in 

In the S-wave case, we encounter the Co function in box diagram, with the analytic 



formula 



ll|: 



Co\pQ,-pQ,0,mQ,mQ] 



r 1 + e 



1 7l' 

- + 2 

e V 



(19) 



The infrared divergence e is canceled by the IR divergence term in self-energy and vertex 
renormalization constants, and the Coulomb singularity term with ^ pole can be absorbed 
into the wave-function by 



a. 



\RsiO)\'il + A— + Ba 



a 



l + A—)il + Bas) + 0{ai). 



(20) 



In the P-wave case, we have to deal with loop-integrals typically as the following ex- 
pression in the box diagram when taking derivative of the relative momentum on the 
denominator of the propagators 

A(/)P 



P{{l-p^^-m?Y{{l+Pif-w?){{l-pi-p2f-m^) 



(21) 



where 2pi is the momentum of the heavy quarkonium and p2 is the momentum of the photon. 
The contribution which is proportional to the term will be omitted when contracted with 
polarization vector. Using the identity A{1) = {A{1) — A{0)) + A{0), we can separate the IR 
divergence into the second term which will be canceled by other diagrams. Three types of 
integrations will appear here, which are given in the Appendix. 

In treating the first term, with the help of the formula I ■ pi = — ((/ — pi)"^ — mQ))/2 
and Dirac decomposition: l^l'^ = g^^Ii + PiPih + {P1P2 + PiP'i)-^^ + ^2^2-^4' ^t)le 
to evaluate most of terms by using LoopTooIs. Note that to get the correct result of the 
integration f dH ((;_p^)2_^2)2((;_)_p^)2_^2) the small number is in the propagators should be 
kept. And we have checked the independence of the final result on e. 

The analytical method is also performed to calculate the one-loop QCD corrections as a 
cross-check for the numerical results. We find the results of the two different methods are 
in agreement. 



IV. NUMERICAL RESULT AND DISCUSSION 

We choose = 10.6 GeV, m^=l.b GeV, mfe=4.7 GeV, a,(2mj=0.26, a,(2mfe)=0.18 
as inputs. As for the charmonium wave-functions at the origin, we choose the results from 
potential model calculations (see the results of the B — T-type potential in Ref. I21I), which 
are listed in Table I. The results of cross sections for e^e~ —>■ 'j* — 7 + A are listed in 
Table II, where ctqed means the QED result and (Tqcd means the corresponding one-loop 
QCD correction. However, if we extract the wave-functions at the origin from the observed 
charmonium decay (e.g., J/ip ^ e^e~ or 77^ 27) widths using theoretical expressions with 



(without) NLO QCD corrections 181] , then the obtained QED cross sections for e^e~ — >■ 
7* — > 7 + A will be larger (smaller) than the S-wave results given in Table II. These are 
the uncertainties due to long-distance matrix elements, and our result in Table II is a rather 
moderate one.*^ 

We see that in most cases the one-loop QCD corrections are negative and moderate, 
except for the Xc2 case, in which the correction is large and is about —75% of the QED result. 
This is very different from the case of double charmonium production e+e^ J/ip + A, 



In Ref.[8| the authors get larger values by using the rjc 27 width with NLO QCD corrections as inputs. 
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TABLE I: Numerical values of the radial wave functions at the origin \R^i 
calculated with the QCD (BT) potential in Ref. 21]. 



for cc and bb 



States 


cc 


bb 


IS 


O.SlGeV^ 


6.477Gel/^ 


2S 


0.529GeV^ 


3.234GeF3 


3S 




2.474GeF3 


IP 


0.075GeV^ 


1.417Gey5 


2P 


0.102Gey^ 




ID 


OMBGeV'^ 





TABLE IL QED results for e+e^ — > 7 + X and the one-loop QCD corrections with wic = 1.5GeV, 
rrih = 4:.7GeV, as{2mc) = 0.26, as{2mh) = 0.18, where oqed means the QED result and (Jqcd 
means the corresponding one-loop QCD correction. 



process 






It 

Vc 




Vb 




II 

Vb 


(yQED{^) 


59.1 


38.6 


33.2 


1.08 


2.19 


0.16 


0.01 


o-QCD (fb) 


-12.5 


-8.19 


-7.04 




-0.55 


~ 


~ 



process 


XcO 


Xcl 


Xc2 


XcO 


xli 


X'c2 


Xho 


Xbi 


Xb2 


crQED{ik>) 


1.66 


18.6 


7.35 


2.25 


25.3 


10.0 


0.46 


2.69 


3.55 


o-QCi?(fb) 


0.28 


-5.13 


-5.49 


0.38 


-6.98 


-7.47 


-0.20 


-0.97 


-1.38 



where one-loop QCD corrections are positive and large in most cases (see Refs.(lO| for LO 



and Refs. 
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121 ] for NLO corrections). 
We find that one-loop QCD corrections do not change the angular distributions of Xco and 
?7c, which read (1 + cos^(6')), confirmed by the effective Lagrangian method. However, when 
including one-loop QCD corrections, the angular distributions of Xci and Xc2 are changed 
from (1.38 + cos2(^)) and (2.72 + cos2(^)) to (1.41 + cos2(^)) and (L61 + cos2(^)), which are 
shown in Fig. [3] and Fig. [4], respectively. 

Since the production of bb mesons is near threshold, we make up a factor in the phase space 
(i-(4m^/I)2)3 for Vb production (P-wave process), and by \if/^¥/(^^2 for Xbj production 



-(4m^s)2 



(S-wave dominated process) as a rough remedy to the phase space integrals. The states of 
rj'fj and rj'fl have not been observed yet, so we use the masses of observed T{2S) and T(35') 
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FIG. 4: Angular distributions for Xci + 7 and f/c + 7 productions in e'^e annihilation up to order 

for replacement. Because of the suppression from the small phase space, the cross sections 
for r]'f^, Vjl and Xfeo are negligible and not useful phenomeno logically. We also choose different 
values of as and vrtc as inputs for comparison, and the obtained cross sections of QED with 
one-loop QCD corrections are shown in Table III. 
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TABLE III: Cross sections of QED with one-loop QCD corrections for varying and m^. 





as = 0.26 


as = 0.21 


mc{GeV) 


TTT-c = 1.4 


rric = 1.5 


rn-c = 1.6 


nic = 1.4 


TJlc = 1.5 


rric = 1.6 




51.3 


46.6 


42.5 


53.7 


48.9 


44.9 




33.5 


29.4 


27.7 


35.1 


31.4 


29.3 




28.8 


26.1 


23.9 


30.2 


27.3 


25.2 


(T(Xco)(fb) 

\/V. »--'-'/ \t/ / 


2.55 


1.94 


1.53 


2.48 


1.81 


1.48 




16.6 


13.5 


11.1 


17.7 


14.6 


12.0 


<Xc2){fh) 


2.31 


1.80 


1.65 


3.53 


2.86 


2.55 




3.46 


2.63 


2.07 


3.36 


2.56 


2.01 


<^{Xci)m 


22.6 


18.4 


15.1 


24.1 


19.7 


16.3 


<^{xc2)m 


3.14 


2.53 


2.25 


4.80 


3.97 


3.48 



From our results, we see the production cross sections for rjc, rj'^, and t]" are about 47fb, 
29fb, and 26fb respectively for rric = l.SGeV and Og = 0.26. Since the new state X(3940), 
which is seen in the spectrum recoiling against the J/ip in the inclusive process e^e~ — > 
J /■ip+ any thing by Belle 22|, is widely believed to be the rj" state (see, e.g., 5j] for discussions). 



we expect it could be seen in the recoil spectrum against the photon. We also find that the 
cross sections for Xci and x'ci 13fh and 18fb respectively, which are much larger than 
that produced in the double charmonium process e~^e~ — > J / + XaiXci) recoiling against 
J/ip. As long as the background of ISR can be largely removed, the exclusive process 
e+e^ — > 7 + Xci(Xci) can be an alternative probe to the Xci meson as well as Xci(2P) 
meson. It will be interesting to see whether there will be signals of X(3872) or X(3940) 
as the candidates of Xci(2P). Whereas the predicted production rates of Xco and its radial 
excitations in e~^e^ — > 7 + Xcoix'co) are much smaller than that in e+e^ J jip + Xcoix'co)- 
By comparing the measurements of the J/ ip process with the photon process, we may clarify 
whether the X(4160), which is copiously produced in association with J/ip, is a radially 
excited state of Xco (say Xco(3-P)), or it is the radial excitation of rjc (say 77^(35')) (see Ref. [5( 
for discussions on the X(4160)). 

We see that although the photon and J/ip meson have the same quantum number of 
J^^ = 1 , when P-wave charmonium states are produced in association with the photon 
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or J jil) in e^e~ annihilation, the behaviors of e^e~ J /ip + X and e^e~ ^ 7 + X are very 
different. In the J/tp case, the associated production of Xco state is prominent, whereas the 
photon favors being associated with the Xci state. Hopefully, the measurement of structures 
recoiling against the photon in the e^e~ — > 7 + X process, especially via the exclusive 
channels J/ipn^TT^ ,ip{2S)'j, J/ip'y and J/ipu, will provide a possible way to search for the 
new heavy quarkonium states, when more experimental data are accumulated in the future, 
and the background from ISR process is largely removed. 
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Appendix 

When we evaluate the numerical result, there are some basic loop integrals. They are 
i;iven by 

[ ^4, {^■P2 f ^ . . 

J P{{l-vi?-m^)\{l+p^f-ml){{l~p^-p2?-ml) ^ > 

/ '^'^WTi V2 2M(i^ %^'^\^2(n V2 ^ = -14-^ (22b) 

/ ^'^mn V2 2^fn^ ^2"^' 2V(, V2 ^ = '^-^^^^ (22c) 

J l\{l-piY -^c)W+Pl) - "^c)((i -P2Y - ^c) 
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where we chose s = ,mc = 1.5GeV. 
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